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Abstract 

Background: The control of Mycobacterium tuberculosis (Mtb) infection begins with the recognition of 
mycobacterial structural components by toll like receptors (TLRs) and other pattern recognition receptors. Our 
objective was to determine the influence of TLRs polymorphisms in the susceptibility to develop tuberculosis (TB) 
in Amerindian individuals from a rural area of Oaxaca, Mexico with high TB incidence. 

Methods: We carried out a case-control association community based study, genotyping 12 polymorphisms of 
TLR2, TLR4, TLR6 and TLR9 genes in 90 patients with confirmed pulmonary TB and 90 unrelated exposed but 
asymptomatic household contacts. 

Results: We found a significant increase in the frequency of the allele A of the TLR9 gene polymorphism rs352139 
(A>G) in the group of TB patients (g.f. = 0.522) when compared with controls (g.f. = 0.383), (Pcorr = 0.01, 
OR = 1.75). Under the recessive model (A/G + A/A vs G/G) this polymorphism was also significantly associated with 
TB (Pcorr = 0.01, OR= 2.37). The association of the SNP rs352139 was statistically significant after adjustment by age, 
gender and comorbidities by regression logistic analysis (Dominant model: p value = 0.016, OR = 2.31; Additive 
model: p value = 0.023, OR = 1.68). The haplotype GAA of TLR9 SNPs was also associated with TB susceptibility 
(Pcorr = 0.02). Differences in the genotype or allele frequencies of TLR2, TLR4 and TLR6 polymorphisms between TB 
patients and healthy contacts were not detected. 

Conclusions: Our study suggests that the allele A of the intronic polymorphism rs352139 on TLR9 gene might 
contribute to the risk of developing TB in Mexican Amerindians. 
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Background 

Tuberculosis (TB) is a leading cause of death worldwide 
[1]. World Health Organization (WHO) estimates one- 
third of the world's population is infected with Mycobacter- 
ium tuberculosis (Mtb) [2] and between the 2000 and 2020, 
about one billion people will become infected and 200 
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million people will develop active TB. Only 5 to 10% of the 
infected individuals develop the clinically active disease in 
their lifetime and the other 90% remain as latently Mtb 
infected individuals [2-4]. The progression to active TB is 
the result of the interplay between environmental, host 
genetic factors and pathogenic characteristics of the Mtb 
strain [3-6]. 

Multiple genes have been involved in the control of Mtb 
and progression to TB [7,8]. In this context, toll like recep- 
tors (TLRs) are a family of phylogenetically conserved 
genes, which are essential for recognition of a broad reper- 
toire of pathogen-associated molecular patterns (PAMPs) 
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on macrophages and dendritic cells and play an important 
role in the innate responses against Mtb [9-13]. 

Genetic variations of TLR1, TLR2, TLR4, TLR6 and 
TLR9 have been associated with the susceptibility to TB 
in different ethnic groups [14-20]. In contrast, other 
studies have failed to demonstrate significant associa- 
tions of TLRs polymorphisms with TB [21-24]. 

To our knowledge, no previous studies have addressed 
the prevalence of TLRs polymorphisms in Mexican 
patients with TB. Therefore, we examined whether poly- 
morphisms in TLR2, TLR4, TLR6 and TLR9 are associ- 
ated with the susceptibility to pulmonary TB in Mexican 
individuals from a rural area with high incidence of TB. 

Materials and methods 

Subjects 

Samples from 180 unrelated individuals (90 patients with 
diagnosis of pulmonary TB and 90 household healthy 
controls) were obtained. The diagnosis of pulmonary TB 
was based on the WHO criteria with presence of clinical 
symptoms, detection of acid-fast bacilli in sputum smear 
samples, Mtb positive cultures in Lowenstein- Jensen 
medium and X-ray evidence of cavitary lesions in lung. 
Only patients older than 18 years were included in the 
study. As a control group, we included unrelated indi- 
viduals that were in close contact with TB patients, all of 
them were asymptomatic and no evidence of positive 
Mtb cultures or radiological lesions in lung were 
detected. Both, TB patients and controls were recruited 
from the programs of TB detection and control in the 
State of Oaxaca and were from the Mazatecan ethnic 
group and were living in the Town called Temascal, a 
rural area near from the City of Tuxtepec in Oaxaca 
State (a high pulmonary TB incidence area located in 
the southeast of Mexico). All studied individuals; their 
parents and grandparents were born in the Mazatecan 
area in Oaxaca State, Mexico. The clinical and demo- 
graphic characteristics of TB patients and controls are 
shown in Table 1. 

The Institutional Review Board (IRB) of the National 
Institute of Respiratory Diseases (INER) reviewed and 
approved the protocol under which all subjects were 
recruited. All subjects provided written informed con- 
sent for these studies, and they authorized the storage of 
their DNA samples at INER repositories for this study. 

DNA isolation and TLR 2, 4, 6 and 9 SNP genotyping 

Genomic DNA was isolated from EDTA-anticoagulated 
peripheral blood by using Qiagen blood mini kit (Qiagen, 
Chatsworth, CA). A total of 12 SNPs at the TLR2 (rs38 
04099, rs7656411, rsl21917864, and rs5743708), TLR4 
(rs5030710, rs5030725, rsl2344353 and rs4986790), TLR6 
(rs6815827) and TLR9 (rs352139, rs5743842, and rs5743 
836) genes were analyzed (See details in Table 2). The 



Table 1 Demographic and clinical characteristics of 
patients and controls 





TB Patients 
(N = 90) 


Controls 

(N = 90) 


P value* 


Age 


46.9 ± 17.8 


42.9 ± 15.9 


ns 


Gender 








M (%) 


38 (42.2) 


32 (35.5) 


ns 


F (%) 


52 (57.7) 


58 (64.4) 


ns 


Comorbidities 








Type 2 Diabetes n (%) 


16/90 (17.7) 


0/90 (0) 


<0.05 


Systemic hypertension n (%) 


4/90 (4.5) 


2/90 (2.2) 


ns 


Immunodeficiency n (%) 


1/90 (1.4) 


0/90 (0) 


ns 


Autoimmunity n (%) 


1/90 (1.4) 


0/90 (0) 


ns 



Data are expressed as means ± SD, means and percentage or n 
and percentage. 

^Comparisons of the continuous variables among groups were performed 
using the Mann-Whitney-Wilcoxon test, whereas the X 2 test was used to 
compare frequencies. 



selection of these SNPs was based on their possible func- 
tional effect and due to existence of previous associations 
with infectious diseases. The SNPs genotyping was per- 
formed using the TaqMan 5' nuclease assay. The PCR 
reaction was carried out using 15 ng of DNA, 12.5 ul of 
2x TaqMan Universal PCR Master Mix, 0.625 ul of 40x 
Assay Mix and 8.8 ul of ddH 2 0. The genotyping was 
performed with validated TaqMan assays (Life Technolo- 
gies/Applied Biosystems, Foster City, CA) TLR2 (C 

22274563_10, C_29420880_10, Exp:06/2017 30 F04 and 

C_27860663_10), TLR4 (C_25761144_10, C 270405 

0_10, C_31784003_10 and C_11722238_20), TLR 6 (C_ 

Table 2 Characteristics of the studied SNPs at TLR2, TLR4, 



TLR6 and TLR9 genes 


Gene 


Chr 


SNP ID 


Location 


Genotype 
(Phenotype) 


TLR2 


4 


rs3804099 


Exon (synonymous) 


C> T (Asn199Asn) 






rs765641 1 


3' near gene 


G> T 






rs121917864 


Exon 

(non-synonymous) 


C> T (Arg677Trp) 






rs5743708 


Exon 

(non-synonymous) 


A> G (Arg753Gln) 


TLR4 


9 


rs5030710 


Exon (synonymous) 


C> T (Ser105Ser) 






rs5030725 


Intron 


G> T 






rsl 2344353 


Intron 


C>T 






rs4986790 


Exon 

(non-synonymous) 


A/G (Asp299Gly) 


TLR6 


4 


rs6815827 


3' near gene 


C> T 


TLR9 


3 


rs352139 


Intron 


A> G 






rs5 743842 


Exon 


C> T (Arg5Cys) 






rs5743836 


5' near gene 


C> T 
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30687096_10) and TLR9 (C 2301953_10, C_308101 

01_10 and C_32645383_10). PCR conditions were: Hold 
95°C/10 min, followed by 40 cycles of 92°C/15 s and 60°C/ 
1 min. All PCR's were performed using 96-well plates in 
a Step One plus real time PCR system {Life Technologies/ 
Applied Biosystems, Foster City, CA). 

Statistical analysis 

Demographic and clinical variables between TB patients 
and controls were analyzed with the Stata statistical soft- 
ware v8.0. We used a significance level of p <0.05. Power 
calculation showed that this significance level would 
yield a power of 80% with a sample size of 90 individuals 
per group. Hardy- Weinberg equilibrium was tested for 
all genotypic combinations of each SNP in TB patients 
and controls. The differences in the distribution of the 
allelic, genotypic and haplotype frequencies of TLR2, 
TLR4, TLR6 and TLR9 polymorphisms in TB patients 
and controls were evaluated by the Mantel-Haenszel, 
Chi-square test with 2x2 contingency tables using the 
EPIINFO statistical software v6.04b. p values resulting 
from the association analysis were corrected by the 
Bonferroni method. Relative risks (RR) with 95% confi- 
dence interval (CI) were estimated as the odds ratios 
(OR). The significance of the SNPs associated to TB was 
adjusted by age, gender and type 2 diabetes (T2D) by lo- 
gistic regression analysis and the associations were tested 
under dominant and additive models. Pairwise LD was 
calculated by the genotype correlation coefficient (r 2 ). 
For all pairs of autosomal SNPs, r 2 measures were calcu- 
lated using the software Haploview v4.2 [25]. 

Results 

Demographic and clinical characteristics of TB patients and 
controls are summarized in Table 1. Prevalence of T2D was 
significandy higher in TB patients (17.7%) when compared 
to controls (0%, p < 0.05). Mean age, gender and the preva- 
lence of other comorbidities such as systemic hypertension, 
immunodeficiency and autoimmunity was similar among 
groups. 

The gene and genotype frequencies of the TLR2, TLR4, 
TLR6 and TLR9 polymorphisms are presented in Table 3. 
Significant deviations from the Hardy- Weinberg equilib- 
rium in the distribution of the TLRs SNP genotypes in TB 
patients and controls were not detected with the exception 
of the SNP rsl21917864 of TLR2 (X 2 = 88, p < 0.05 in both 
TB patients and controls). 

We found a significant increase in the frequency of the 
allele A of the TLR9 gene polymorphism rs352139 (A>G) 
in the group of TB patients (gene frequency (g.f.) = 0.522) 
when compared to controls (g.f. = 0.383, Pcorr = 0.01, 
OR = 1.75, 95% CI = 1.15-2.67). Under a dominant model 
(A/A + A/G vs G/G), this polymorphism was associated 
withTB (Pcorr = 0.009, OR = 2.37, 95% CI = 1.22-4.60). 



Table 3 Distribution of TLR2, TLR4, TLR6 and TLR9 gene 
polymorphisms in patients with TB and controls 

TB patients Controls 
(N=90) (N=90) 
SNP Genotype n F n F Pcorr 



Allele 


rs3804099 


T/T 


59 


(0.656) 


48 


(0.533) 


ns 




C/T 


26 


(0.289) 


36 


(0.400) 


ns 




C/C 


5 


(0.056) 


6 


(0.067) 


ns 




T 


144 


(0.800) 


132 


(0.733) 


ns 




C 


36 


(0.200) 


48 


(0.267) 


ns 


rs765641 1 


T/T 


58 


(0.644) 


53 


(0.589) 


ns 




G/T 


29 


(0.322) 


33 


(0.367) 


ns 




G/G 


3 


(0.033) 


4 


(0.044) 


ns 




T 


145 


(0.806) 


139 


(0.772) 


ns 




G 


35 


(0.194) 


41 


(0.228) 


ns 


rs121917864 


C/C 


- 


- 


- 


- 






C/T 


90 


(1.0) 


90 


(1.0) 


ns 




T/T 


" 


" 


" 


" 






C 


90 


(0.500) 


90 


(0.500) 


ns 




T 


90 


(0.500) 


90 


(0.500) 


ns 


rs5743708 


G/G 


90 


(TO) 


90 


(TO) 


ns 




G/A 


- 




- 








A/A 


- 




- 








G 


180 


(1.0) 


180 


(1.0) 


ns 




A 


■ 


- 


- 


- 




rs5030710 


T/T 


87 


(0.967) 


87 


(0.967) 


ns 




C/T 


3 


(0.033) 


3 


(0.033) 


ns 




C/C 


- 


- 


- 


- 






T 


1// 


(0.983) 


1// 


(0.983) 


ns 




C 


3 


(0.017) 


3 


(0.017) 


ns 


rs5030725 


T/T 


86 


(0.978) 


89 


(0.989) 


ns 




T/G 


2 


(0.022) 


1 


(0.011) 


ns 




G/G 


- 


- 


- 


- 






T 


178 


(0.989) 


179 


(0.994) 


ns 




G 


2 


(0.011) 


1 


(0.006) 


ns 


rs 12344353 


T/T 


89 


(0.989) 


90 


(1.0) 


ns 




C/T 


1 


(0.011) 






ns 




C/C 














T 


179 


(0.994) 


180 


(1.0) 


ns 




C 


1 


(0.006) 






ns 


rs4986790 


A/A 


88 


(0.978) 


89 


(0.989) 


ns 




A/G 


2 


(0.022) 


1 


(0.011) 


ns 




G/G 














A 


178 


(0.989) 


179 


(0.994) 


ns 




G 


2 


(0.011) 


1 


(0.006) 


ns 
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Table 3 Distribution of TLR2, TLR4, TLR6 and TLR9 gene 
polymorphisms in patients with TB and controls 

(Continued) 

TLR6 rs6815827 



TLR9 rs352139 



rs5743842 



rs5743836 



c/c 


88 


(0.978) 


88 


(0.989) 


ns 


C/T 


2 


(0.022) 


1 


(0.011) 


ns 


T/T 


- 


- 


- 


- 




C 


178 


(0.989) 


1// 


(0.994) 


ns 


T 


2 


(0.011) 


1 


(0.006) 


ns 


G/G 


19 


(0.211) 


35 


(0.389) 


0.01 


A/G 


48 


(0.533) 


41 


(0.456) 


ns 


A/A 


23 


(0.256) 


1-1 


(0.156) 


ns 


G 


86 


(0.478) 


111 


(0.617) 


0.01* 


A 


94 


(0.522) 


69 


(0.383) 


0.01" 


T/T 


_ 


_ 


_ 


_ 




T/C 


2 


(0.022) 






ns 


c/c 


88 


(0.978) 


90 


(1.0) 


ns 


T 


2 


(0.011) 






ns 


C 


178 


(0.989) 


180 


(1.0) 


ns 


T/T 


82 


(0.911) 


78 


(0.867) 


ns 


T/C 


8 


(0.089) 


12 


(0.133) 


ns 


ac 












T 


172 


(0.956) 


168 


(0.933) 


ns 


C 


8 


(0.044) 


12 


(0.067) 


ns 



* OR = 0.42, 95% CI = 0.21-0.81. 
** OR = 0.56, 95% CI = 0.38-0.86. 
*** OR = 1.75, 95% CI = 1.15-2.67. 

ns: non significant; F: frequency; Pcorr: p values were corrected by Bonferroni. 
Statistically significant values are shown in bold. 



The frequency of the homozygous genotype G/G of 
the SNP rs352139 was considerably lower in the group 
of TB patients (g.f. = 0.211) compared to controls (g.f = 
0.389) (Pcorr = 0.01, OR = 0.42, 95% CI = 0.21-0.81). 
Consequently the frequency of the allele G was signifi- 
cantly higher in the control group (g.f. = 0.617) than in 
TB patients (g.f. = 0.478) (Pcorr = 0.01, OR = 0.56, 95% 
CI = 0.38-0.86). We also analyzed the distribution of 
allele and genotype frequencies of the rs352139 poly- 
morphism in males and females from the TB and control 



groups. A significant association of the allele A under a 
dominant model (A/ A + A/G vs G/G) was found with 
TB in females (Pcorr = 0.00008, OR = 5.2 (95% CI = 2.3- 
12.0). 

Importantly, the association of the genotypes of the 
SNP rs352139 with TB was confirmed by logistic regres- 
sion analysis under the dominant (A/A + G/A vs G/G; 
p = 0.01, OR= 2.31, 95% CI= 1.50-4.63) and additive 
(2A/A + G/A, to G/G, p = 0.02, OR = 1.68, 95% CI = 
1.07-2.65) models after adjustment for age, gender and 
T2D, Table 4. 

We also compared the distribution of all TLRs SNPs 
in TB patients without T2D and controls. The allele A 
of the SNP rs352139 was associated to TB in patients 
without T2D (Pcorr = 0.02, OR = 1.63, 95% CI = 1.06- 
2.57). Also, the dominant model (A/ A + A/G vs G/G) 
was associated with TB, excluding individuals with T2D 
(Pcorr = 0.02, OR = 2.26, 95% CI = 1.12-4.55). Signifi- 
cant differences in the distribution of genotypes or 
alleles of the TLR9 polymorphisms (rs5743842 and 
rs5743836) were not detected. 

The allele and genotype frequencies of the TLR2 
(rs3804099, rs7656411, rsl21917864, rs5743708); TLR4 
(rs5030710, rs5030725, rsl2344353, rs4986790) and TLR6 
(rs68 15827) polymorphisms were similar amongst TB 
patients and controls. 

The distribution of TLR2, TLR4 and TLR9 haplotypes 
in TB patients and controls are listed in Table 5. We 
found a significant association of the haplotype GAA of 
TLR9 with TB (TB patients: hf = 0.522 versus Controls: 
hf = 0.383; Pcorr = 0.02). No significant differences were 
detected in the distribution of TLR2 and TLR4 haplo- 
types between TB patients and controls. In general, low 
LD (r 2 ) values were observed for all TLR2, TLR4 and 
TLR9 polymorphisms, Figure 1. 

Discussion 

Several studies have demonstrated that genetic factors 
play a major role in the susceptibility to TB [7,8]. Evi- 
dence suggests that polymorphisms within TLRs may 
cause an impairment of the immune response against 
TB; possibly modifying the TLRs affinity to Mtb derived 



Table 4 Association of the TLR9 polymorphism rs352139 (A>G) genotypes with pulmonary TB after logistic 
regression analysis 







Genotype frequency (%) 












G/G 


G/A 


A/A 


MAF (A) 


Model 


OR (95% CI) 


P value 


Control (n = 90) 


0.38! 


3 0.456 


0.156 


0.383 








TB (n = 90) 


0.2 V 


1 0.533 


0.256 


0.522 


Dominant (A/A + G/A vs G/G) 
Additive (2A/A + G/A, to G/G) 


2.31 (1.50-4.63) 
1.68 (1.07-2.65) 


0.016 
0.023 



Associations were tested using logistic regression adjusting for Age, Gender and T2D. 

TB: Tuberculosis; MAF: Minor allele frequency. Statistically significant values are shown in bold. 
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Table 5 Haplotype distribution of TLR2, TLR4 and TLR9 
polymorphisms in patients with TB and controls 

Haplotypes* TB patients Controls X 2 Pcorr 



(N = 90) (N = 90) 





hf 


hf 






TLR2 










TTGT 


0.349 


0.324 


0.243 


ns 


TCGT 


0.349 


0.324 


0.243 


ns 


CTGG 


0.046 


0.071 


1.053 


ns 


GCGG 


0.046 


0.071 


1.053 


ns 


CTGT 


0.054 


0.062 


0.103 


ns 


GCGT 


0.054 


0.062 


0.103 


ns 


TTGG 


0.051 


0.042 


0.154 


ns 


TCGG 


0.051 


0.042 


0.154 


ns 


TLR4 










TTTA 


0.972 


0.972 


0.0 


ns 


TTCA 


0.011 


0.017 


0.203 


ns 


TLR9 










GGA 


0.433 


0.550 


4.901 


ns 


GAA 


0.522 


0.383 


7.007 


0.02 


GGG 


0.033 


0.067 


2.105 


ns 



Haplotypes: TLR2 [rs3804099, rs121917864, rs5743708, rs7656411]. 
TLR4 [rs12344353, rs5030725, rs5030710, rs4986790], 
TLR9 [rs5743842, rs352139, rs5743836]. 

ns: non significant, hf: haplotype frequency. Statistically significant values are 
shown in bold. 



ligands leading to alterations in the signal transduction 
of anti-Mtb immune responses [13,16]. These studies 
have been performed in different populations including 
Asians, African-Americans and Caucasians, but the 
results of genetic association studies are controversial 



[14-24]. Recent studies suggest that TLRs polymor- 
phisms might not contribute markedly to TB susceptibil- 
ity [21,26]. However, the current study provides evidence 
that TLRs polymorphisms contribute to pulmonary TB 
susceptibility in Mexicans. We found a significant associ- 
ation of the allele A of the SNP rs352139 (A>G) located 
in the intronic region of TLR9, on chromosome 3, with 
the susceptibility to TB in patients from a rural area in 
the State of Oaxaca in the southeast of Mexico. The asso- 
ciation of this SNP was confirmed by logistic regression 
analysis, adjusting for age, gender and T2D. A significant 
association with the dominant A/ A + A/G vs GG with 
TB was also detected and particularly with female. A de- 
creased frequency of the homozygous genotype G/G of 
the SNP rs352139 was detected in the group of TB 
patients, suggesting a protective role of the allele G. 
Moreover, the TLR9 haplotype containing the allele A 
of this SNP (Haplotype GAA) was associated with the 
susceptibility to TB. 

TLR9 is expressed by different cells of the immune 
system and is essential in the activation of the innate 
responses against Mtb infection [9]. TLR9 recognizes 
unmethylated CpG DNA motifs and functional studies 
revealed that the binding of TLR9 is necessary to drive 
the Thl immune response [27]. The importance of 
TLR9 in the antimycobacterial responses has been sup- 
ported in experimental models using the TLR9-deficient 
mice [9,10]. 

In line with our results, the allele A of the SNP 
rs352139 has been strongly associated with susceptibility 
to TB in African-Americans but not with Caucasians 
and African patients from Guinea-Bissau [17]. Func- 
tional studies suggest that the intronic TLR9 polymorph- 
ism rs352139 A>G exerts a regulatory effect of TLR9 
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gene expression. The presence of G allele and the C 
allele of the promoter polymorphism -1486, is associ- 
ated with downregulation the TLR9 expression at the 
transcriptional level. In contrast, the presence of the 
allele A (in presence of the allele T of the promoter poly- 
morphism T allele at -1486) induces the up-regulation of 
TLR9 expression [28]. It is possible that this variation influ- 
ences signaling by creating an alternative splicing site, 
affecting the TLR9 mRNA transcription. Otherwise, the 
rs352139 SNP could be a likely marker in LD with a poly- 
morphic regulatory region that controls TLR9 expression 
or a functional coding region polymorphism [28]. However, 
the effect of this intronic SNP in the induction or down- 
regulation of TLR9 expression is not well defined, particu- 
larly in the context of infectious diseases. 

Different from other pathogens, Mtb infection can 
persist in the host for long periods in a dormant or latent 
state, even in a fully functioning immune system [29]. TB 
immunity is mediated by Thl-type responses; nevertheless 
this response is apparendy not enough to definitively 
eradicate the chronic infection or the latent bacilli [30]. 
It is possible that the presence of the allele A of the 
rs352139 SNP in individuals with latent TB infection may 
result in a high TLR9 expression that contribute to the in- 
duction of high levels of proinflammatory cytokines. A 
pro-inflammatory environment in lung confers protection 
against Mtb infection. However, a pathogenic or exacer- 
bated pulmonary inflammation may exert deleterious ef- 
fects [31], promoting the infiltration inflammatory cells to 
the lung that become Mtb infected. A dysregulated inflam- 
matory environment, potentially can also contribute to the 
tissue damage and dissemination of the mycobacteria to 
other areas of the lung and peripheral tissues. It is import- 
ant to remark that the allele A of this polymorphism of 
TLR9 has been associated with the susceptibility and se- 
verity of parasitic and viral infections [32,33]. 

Interestingly, the polymorphism rs352139, located in 
the intronic region of TLR9, has been associated with 
the susceptibility to TB in Indonesian females [14]. 
Besides, it is noteworthy that this polymorphism was 
associated with pulmonary TB in females from South 
India [34]. In our study, we found a significant associ- 
ation of this polymorphism, particularly the homozygous 
A/ A and heterozygous A/G genotypes, with the suscepti- 
bility to TB in females. The reason why these genotypes 
contribute to the susceptibility to TB in females is not 
fully clear. In this regard, previous studies have sug- 
gested that significant differences exist in TB incidence 
between males and females [35,36]. The physiological 
hypothesis of infectious diseases susceptibility, suggest 
that differences in the expression of sex hormones and 
genetic variability between males and females might 
contribute in the differential incidence of infectious dis- 
eases among them [37]. Moreover, in vitro studies have 



revealed that progesterone regulates the production of 
type I interferon (IFN-a), possibly through the interference 
of MyD88-dependent activation of interferon regulatory 
factor 7 (IRF-7) after stimulation with the TLR9 ligand 
CpG in both, mice and humans. Consequently, the inhibi- 
tory effect of this female sex steroid hormone on IFN-a 
production may impair inflammatory responses and pat- 
tern recognition receptor signaling against pathogens [38]. 

From the population genetics perspective, it is interest- 
ing that this SNP is associated with susceptibility to TB 
across Asian and Amerindian groups. The Mazatecans 
inhabit the mountains of northern Oaxaca (Mexico). 
They belong to the ancient Central-American cultural 
groups and may be related to the Olmecs (the most 
ancient Meso-American culture). The first Amerindian 
Natives are believed to have come from Asia through 
the Bering land bridge between 30,000-12,000 years 
before the present. These assumptions have been based on 
cultural, morphological and genetic connections between 
Native American and Asian populations [39-41]. Sup- 
porting the hypothesis of the presence of Asian genes in 
Amerindian populations, in previous published studies 
we have detected Asian haplotypes (HLA-DRB1*15:02- 
DQB1*0602) in Mazatecans [42]. Also, a short genetic 
distance between Chinese and Native Americans was 
detected using Alu insertions data [43]. Anthropological 
genetic studies have suggested that some Asian genes may 
have been introduced to Native Americans by a Trans- 
Pacific route of migration from South East Asian popula- 
tions [44,45]. 

Our study has some limitations, including its relatively 
small sample size, restricted by the study's focus on 
patients from Amerindian ancestry with a recent diagno- 
sis of pulmonary TB. In addition, we chose not to in- 
clude individuals with familial history of autoimmune 
and neoplastic disorders. The second limitation is the 
lack of functional assays to determine the effect of TLR9 
variation in the susceptibility to pulmonary TB. The 
strengths of our study includes the analysis of the influ- 
ence of TLRs polymorphisms in the susceptibility to TB 
in a well genetically characterized ethnic group of Mexi- 
can Amerindians and the replication of the associations 
of the TLR9 intronic polymorphism previously described 
in other Asian ethnic groups. 

Conclusions 

In summary, our findings suggest that the allele A of the 
polymorphism rs352139, located in the intronic region of 
the TLR9 gene could be implicated in the susceptibility to 
TB in Mazatecan Amerindians from the Oaxaca State in 
Mexico. Additional studies in a large number of TB pa- 
tients and functional studies may help to establish the true 
significance and the possible deleterious effect of this poly- 
morphism in the susceptibility to pulmonary TB. 
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